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Abstract 



The evolution of the EMC effect as a function of atomic mass A is considered for the first 
time for the lightest nuclei, D, 3 H, 3 He and 4 He, with an approach based on the Bethe-Salpeter 
formalism. We show that the pattern of the oscillation of the ratio r A (x) = F^/F^ varies with 
A, unlike the pattern for nuclei with masses A > 4, where only the amplitude of the oscillation 
changes. 
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The understanding that nucleon structure can not be regarded as unrelated to nuclear struc- 
ture has been the main outcome of the European Muon Collaboration (EMC) Q. The nuclear 
environment modifies the nucleon partonic structure in such a way that the ratio of the nuclear 
and deuteron structure functions, r (x)= F A /F^ , deviates from unity, resembling an oscillation 
^ with respect to the line r A (x) = 1. Even though the A dependence of the ratio r A {x) could be 

explained with conventional nuclear structure considerations in the range A > 4 £|, the origin of 



Qh! the effect remained obscure ||. 

The argument that the origin of the EMC effect is closely related to saturation of the short- 



range binding forces in 3- or 4- nucleon systems was first presented in j|, || . Previous attempts to 
exploit this physics for explanation of the effect failed because of both theoretical difficulties and 
a lack of understanding on how the saturation would show up in the observables. This excuses 
the statement of ||, that the data on r (x) do not directly correlate with the binding energy 
per nucleon. The saturation, according to ||, had to manifest itself not in the amplitude of the 
oscillations, but in the pattern of the x dependence of r A (x), namely in the positions of the three 
cross-over points Xi, i = 1 — 3, in which r A (x) = 1. Such a pattern can be clearly seen from new 
data of SLAC || and NMC . There does not exist any data on the EMC effect in the range of 
A < 4. Most challenging therefore is to evaluate how the effect evolves in the range of the lightest 
nuclear masses, r D (x) — ► r A=3 (x) — > r A=4 (x), starting from the basic property of nucleons to 
form bound states. 

In the present Letter, we perform derivations of the relative changes in the nuclear structure 
function F A (x) with respect to the isoscalar nucleon one, F^(x) = ^(F^x) +i ? 2 1 (x)), where p and 
n denote the free proton and the free neutron respectively. On the other hand, the comparison 
with experimental data can be done only in terms of r (x), obtained with the deuteron structure 
function F^{x) as a reference. Therefore our final results will be presented for both cases. In the 
considered range of x ( 0.3 < x < 0.9 ) the experiments (see Ref. [||]) are consistent with no Q 2 
dependence of r A (x). We perform numerical calculations for a fixed Q 2 of 10 GeV 2 . 
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Our approach originates from the Bethe-Salpeter formalism j8), which allows one to treat 
nuclear binding effects by using general properties of nucleon Green functions B. The approach 
is model independent in the sense that it does not require any assumption about the nuclear 
structure except that the nuclear fragments have a small relative energy. This allows us, as has 
been shown previously |J, to derive the deuteron structure function in the form: 

r,n/ x f d3k ( nN / s/ h\ M D -2E (fF 2 N (x N ) \ T 9/ , x 

where E is the on- mass-shell nucleon energy E 2 = k 2 + m 2 , and k is the relative three momentum 
of the bound nucleons, is the mass of a deuteron. The nucleon Bjorken variable is defined 
as in = x-£,m/(E — k$). The function ^f 2 (k) is an analog of the three dimensional momentum 
distribution and is directly related to the Bethe-Salpeter (BS) vertex function for the deuteron 

r D (p,k) = {r°(p,% o ^_^ : 

* 2(k) = iM^ib 1 ^ • (2) 

The two-term structure of Eq. ([l]) defines a deviation of the ratio F 2 ) / P<^ from unity, which 
is generally considered as a signature of the modification of the nucleon structure. As has been 
shown in |J, the second term results from the relative time dependence of the amplitude of 
lepton deep inelastic scattering (DIS) off a bound nucleon. By extending the approach for light 
nuclei, A = 3, 4, we have discovered that it is the relationship between similar terms, which is 
responsible for evolution of the nucleon partonic structure. 

The hadronic part of the DIS amplitude (hadronic tensor) is related with the forward Compton 
scattering amplitude T^ u by using the unitarity condition 



W£(P,q) = ±ImT£,(P,q), (3) 

and T^ v is defined as a product of electromagnetic currents averaged over nuclear states, 

l£(P, q) = ij d A xe^(A\T (J, (x) J u (0)) \A). (4) 

Starting from a field theory framework one can define the matrix element in terms of solutions of 
the n-nucleon Bethe-Salpeter equation and n-nucleon Green functions with the insertion of the 
T-product of electromagnetic currents: 

T*(P, q)= J dKdK'T(P, K,)S n {P K)G 2(n+l) ^{q; P, K, K.')S n (P, K')T(P, £'), (5) 

where K, denotes a set of momenta, which describes relative motion of nucleons, K. = ki, . . . , k n —±, 
dJC = d 4 /c n /(2-7r) 4 . . . d 4 fc n _i/(2-7r) 4 , and P is the total momentum of the nucleus. The BS vertex 
function T(P,IC), introduced to describe a nucleus in terms of virtual nucleon states, satisfies the 
homogeneous Bethe-Salpeter equation 

T(P,K) = -[ dK'G 2n (P,lC,lC')S n (P,}C')T(P,lC'). (6) 

The S n (P,K') is direct product of n-nucleon propagators. The Cr2 n (P, /C, K!) term denotes the 
irreducible truncated n-nucleon Green function which is defined as follows: 

G 2ri ,(P,/C,/C') = S n -\P,1C)5{1C - K,') - G 2n {P,KX), (7) 
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where G 2n is an exact n-nucleon Green function. The Green function (?2(n+l) „ represents Comp- 
ton scattering of a virtual photon off a system of re-virtual nucleons. All irreducible interaction 
corrections to the imaginary part of T^ v are suppressed by powers of 1/(Q 2 ) §§. This justifies 
consideration of the zeroth order term of G 2 (n+x) v ' 

G 2{n+l)ixv {q; P, K) = £ O^fa; P, k^S^h, . . . k^, k i+1 , . . . fc n _ 1 )<5(/C-X: / )+0(l/Q 2 ), (8) 

i 

where G^ v is connected with nucleon Compton amplitude: G4 /xu (q; P, k{) = u{k.\)T^ v {q, fej)u(k;). 
Then T^ u can be rewritten in terms of the off-mass-shell T^: 

T^(P,q) = J dlC^T^qMk^SniP^iMk^Ti^^Sni^lC^i^JC), (9) 

* 

Integration over k{Q can in principle relate T^ v with on-mass-shell nucleon Compton amplitude. 
This can be realized only after the singularities in nucleon propagators and the BS vertex functions 
are taken into account || . Unlike the deuteron case, where singularities in the BS vertex function 
can be neglected, in this case they are connected with nucleon-nucleon bound states, which lie 
in the range of low relative momenta. One can express the singularities explicitly by introducing 
the "bare" BS vertex function 7, which is regular with respect to the relative nucleon momenta: 

r(p,/c) = - J d)Cg2n{P,JC,)c')s n (p,K:'MP,}c'), (10) 

where c/ 2n denotes the regular part of n-nucleon Green function at P 2 — > M|. This function, 
however, contains singularities of m-nucleon Green functions (m < re). For example, in case of 
3 He the function g§ contains a pole of G4, which corresponds to a bound deuteron and nucleon- 
nucleon continuous spectrum (74: 



G 4 [ — + k,ki,k , 1 ) = ^- 2 3 —+54 ( — + k,k 1 ,k[) . (11) 

3 ' f ^+k] -M 2 ^ 3 



2P , , , rJ \ T D ( 2 f + k,k 1 )T D (^ + k,k' 1 ) , _ (2P 

f + kf-Ml 

For 4 He one has additionally the 3 He and 3 H poles. Substituting these expressions into Eq. (0) 
and integrating over the relative energy of different nuclear fragments we derive the 3 He, 3 H and 
4 He Compton amplitudes respectively in terms of physical amplitudes of the fragments. Using 
the projection operator g nv and Eq. (Q) one gets 

lim g^W^ A \P,q) = --F 2 N(A) (x) . 

Introducing now Bjorken variables xa = 2 p A . an d = 2 p N - > we ^ n< ^ ^2 f° r H e an d H in the 
form: 

^2 He (^He) = (12) 



d 3 k 



£ p ^p) + ^ *2 + ^ *P ^ p dx D 



^2 H (^H) = ^2 He 0£3He)lp~n 



^3 He (k), 
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and for 4 He in the form: 



(13) 



d 3 k 



+ 



(2tt) 3 

E n 



E„ 



L *?(*p) + 



E: 



3 H 



^(-h) + ^ 



dE£(x p ) 



a 4Hc 



EL 



dx 



E: 



3 He 



Es u , 



^F^(x 3Ue ) + 



a 4Hc 



p -^ 3 H 
4 He 



dF* H (s3 H ) 
X3 H - 



dx 



dFj(x n ) An 



3 H 

dF 2 3He (^ He )' 



3 He 



3 He 



^He(k), 



where Aj^ = —Ma + -En + -EU-i is the binding energy of the corresponding nuclear fragment. 
The three dimensional momentum distributions $^(k) are defined via the "bare" Bethe-Salpeter 
vertex functions. For example for 3 He one has: 



* 3 He(k) 



AE p E D M mc (M D -E p - E D y 



f d*h d% s n J2P 
J (2vr)4(27r) 47 ^ k ' k ^ 2 { 3 + /c ^ 



(14) 



x r 



D 



2P 



+ k,h r 



D (IP 



2P 



E<(k)#) 7 3He 0P,Mi) 



where fcop = — I s — -Ep- Since presently there are no realistic solutions of the Bethe-Salpeter 
equation for a bound system of three or more nucleons, one has to use phenomenological mo- 
mentum distributions for numerical evaluations. It is reasonable to assume that the momentum 



distributions in Eqs. (12) and ( |i3|) can be related with the distributions extracted from experi- 
mental data. In the calculations we make use of the distributions available from [|n]] and |11[| . 
The contribution arising from continuous spectra (ppn for 3 He and ppnn for 4 He) is small in the 
considered kinematical range and does not change comparison of the final result with the data. 
This justified some simplifications which resulted in rather transparent form of Eqs. (^) and (|l~3|). 
The procedure has been consistently taken into account in the normalization of the momentum 
distributions $3 H and $4 H • 



This result reduces to the one obtained within the x-rescaling model [12] and for A = 3 
becomes: 



X3 He 



/ 



D/ 3 He 



(15) 



M 3Hc 



where e = A p He has the meaning of a nucleon (deuteron) separation energy and / p(D)/3Hc (y,e) 
are the 3 He spectral functions for a bound proton (deuteron): 



/ 



p(D)/ 3 He 



d 3 k » . m .( E p r D \ - k 3 

- yS[y- 



(27T) 



r* 3 He(k) 



E 



P(D) 



S(e- (E P + E D -M* He )) 



Input structure functions F% ( x ) are introduced via parameterizations based on the measure- 
ments of the proton and the deuteron structure functions by DIS experiments. We used the most 
recent parameterization of F^x^Q 2 ) found in |l3j and fixed the value of Q 2 to 10 GeV 2 . The 
structure function F^{x) is evaluated from F^ix) and from the ratio F^ix^/F^ix) determined 
in [p^ |. We have verified that the uncertainties in F% (x) are suppressed in the obtained ratio 
r (x) and thus can be neglected in the considered kinematic range. 

The results of the numerical calculations, which show how the free nucleon structure function 
F^{x) (A = 1) evolves to the deuteron (A = 2) and helium (A = 3 and 4) structure functions, are 
presented in Fig. 1(a). The evolution, which starts from F^fa), is shown in Fig. 1(b). Contrary 
to what is observed for nuclei with masses A > 4, the pattern of the oscillation of r A {x) changes 
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its shape in the range of A < 4. The rate at which the changes occur is consistent with the fast 
buildup of the short range binding forces. 

We compare our results for the ratio F 2 ne (x) / F 2 D (x) with the available data from || in 
Fig. 2. The position of the cross-over point, obtained from our calculations as X3 = 0.919, is in 
reasonable agreement with the extrapolated data. On the other hand, the corresponding point 
for A = 3, x% = 0.855, falls within the interval 0.84 - 0.86, which is where the ratios r A (x) (A > 
4) cross the line r A = 1 (cf. Ref. ||). This means that the pattern of the EMC effect observed 
in such dense nuclei as metals is being reached at A = 3. The larger value of x% at A = 4 must 
be related with the anomalous binding energy of 4 He. 

A fundamental relation follows from the obtained results. Since binding corrections have the 
same form in Eqs. (12) we can write 



1= J^(F 3 2 ^(x)-F^(x)) = J^(FZ(x)-F?(x)). (16) 


The result represents the Gotfried sum /, which has been often studied experimentally from the 
combination of F 2 (x) and F 2 ) (x)(cf. Ref. [|l5|]). Such a combination is equal to I to within a 
correction proportional to F 2 i (x = 0). Indeed, as follows from Eq. (|l]), 

= (2F$(x) - 2F?(x)) = / - 2 < Mp ~J^ d)d if(x = 0) . 



Apparently, such tests cannot be performed rigorously because F 2 (x) is unknown at x = 0. On 
the other hand, an experiment, which used 3 He and 3 H targets, would be able to measure the 
nucleon isospin asymmetry independently of the model uncertainties in the binding corrections. 

In conclusion, the method for the model- free calculations of the evolution of the nucleon 
structure in the lightest nuclei has been developed as the extension of an approach based on 
the Bethe-Salpeter formalism. The method allows one to express F 2 (x) in terms of structure 
functions of nuclear fragments and three dimensional momentum distributions. As a result, 
F 2 A /F 2 ^ have been evaluated numerically without finding solutions of Eq. (0). 

The obtained pattern of distortions of the nucleon structure function prove that the EMC 
effect in the lightest nuclei, D, 3 H, 3 He and 4 He, is basically the manifestation of the short range 
binding forces in the nucleon parton distributions. The quantitative predictions for 3 He and 4 He 
nuclei, which have to be verified in future experiments at HERA or CEBAF, indicate that the 
EMC effect in heavy nuclei can be naturally understood as distortions of the nucleon parton 
distributions in 3 He or 3 H, which are modified by the nuclear density effects. 

Finally, the obtained results prove that in the EMC effect range (0.3 < x < 0.9) the two- 
nucleon interactions can be considered as the dominant mechanism for the description of the 
nuclear binding forces. 

We thank S. Akulinichev, A. Antonov, A.M. Baldin, S. Kulagin and V. Nikolaev for useful 
discussions. A.M. acknowledges the warm hospitality of the Special Research Center for the 
Subatomic Structure of Matter, Adelaide, Australia. This work was supported in part by the 
RFBR grant N96-15-96423. 
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Figure 1: (a) The ratio F£ (x) / F 2 N (x) . (b) The ratio of F^(x) to the deuteron structure function 
FrP(x) (A = 4) and to the combination of structure functions (2F 2 D (x) + F 2 p (a;))/3 (A = 3). The 
dashed curve in FIG. (a) shows the result of calculations, described in the text, for A = 2. The results 
for A = 3, 4 are shown with the solid curves. 
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Figure 2: The ratio F^i^x) / ' F^(x). Results of the calculation, described in the text, are shown with 
the solid curve. The data are from Ref. [5] (filled circles) and Ref. [6] (empty circles). 
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